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ABSTRACT 


Transport  property  measurements  with  ultrasonic  techniques 
are  discussed  over  the  entire  temperature  (100-15,  000  K)  and  pres¬ 
sure  (1-1000  atm)  ranges.  A  chart  of  various  techniques  for  the 
determination  of  transport  properties  in  the  various  temperature  and 
pressure  ranges  is  given.  Sound  absorption  and  sound  speed  deter¬ 
minations  were  used  to  obtain  thermo  physical  properties  of  argon, 
helium,  nitrogen  and  air  at  room  temperature  and  pressures  between 
1-500  atm.  Several  variations  of  pulsed  ultrasonic  techniques  were 
used  to  obtain  maximum  precision  in  the  gas  property  measurements. 
Sound  speeds  accurate  to  approximately  0,  1%  were  obtained  in  all  the 
gases.  Thermal  conductivity  and  specific  heat  ratios  were  deduced 
for  air. 


A  cavity  technique  was  developed  for  high  accuracy  (<  0.  1%) 
viscosity,  thermal  conductivity  and  bulk  viscosity  determinations. 
Initial  tests  were  run  to  proof  the  technique  at  room  temperature  and 
pressure.  Sound  speed,  Prandtl  number  and  viscosity  measurements 
were  made  which  agree  with  accepted  values  to  within  1%. 
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SECTION  I 
INTRODUCTION 


The  objective  of  this  work  is  the  experimental  determination  of 
gas  properties  at  high  temperatures  and/or  pressures,  with  emphasis 
on  transport  properties.  The  immediate  goals  of  the  present  work  are 
to  extend  existing  tabulations  for  immediate  engineering  needs  and  to 
provide  long  range  support  for  theoretical  work.  The  tem'perature  range 
is  300  K  to  20,  000  K,  and  the  pressure  range- is  1  to  10^  atmospheres. 

The  principal  transport  properties  of  interest  are  viscosity,  thermal 
conductivity  and  diffusion.  In  addition,  consideration  of  other  properties 
such  as  electrical  conductivity,  optical  depth  and  certain  thermodynamic 
properties  are  of  interest. 

An  updated  chart  of  the  useful  techniques  in  the  temperature- 
pressure  regions  of  interest,  the  high  pressure  region  (300  <  t <  1000  K, 
p  <  10,  000  atm),  the  chemical  reaction  region  (1000  <  t  <  8000°K, 
p  <  300  atm),  and  the  plasma  region  (800°K  >  t),  as  discussed  in  a  pre¬ 
vious  report  (Ref.  1),  is  presented  in  Fig.  1.  Extensive  ultrasonic  mea¬ 
surements  were  made  at  Panametrics  in  the  plasma  region  (Refs.  2,3,4). 

A  paper  is  in  preparation  in  which  these  measurements  in  nitrogen, 
oxygen  and  air  are  used  to  obtain  transport  properties.  While  there 
have  been  considerable  advances  in  the  theory  of  transport  properties 
of  plasmas  in  recent  years  (Refs.  5,  6),  the  development  of  new  experi¬ 
mental  techniques  seems  to  be  relatively  stagnant.  Some  new  develop¬ 
ments  with  pressure  and  magnetic  probes  are  in  progress  (Ref.  7).  In 
addition,  work  is  being  carried  out  to  extend  transport  property  mea¬ 
surement  in  plasmas  up  to  100  atm  (Ref.  8).  Measurements  of  transport 
properties  of  plasmas  are  of  interest  for  engineering  of  arc  tunnel  test 
facilities,  magneto-hydrodynamic  generators,  and  light  source  development. 

The  ultrasonic  and  Wollaston  prism  interferometric  techniques  in 
the  shock  tube  were  used  at  Panametrics  in  the  chemical  reaction  region. 
The  arc  techniques  do  supply  information  in  the  temperature  range  where 
dissociation  occurs.  However,  the  internal  inconsistency  of  thermal 
conductivity  measurements  from  arcs  is  anomalously  large  in  regions 
of  dissociation  (Refs.  9,  10,  11).  The  authors  feel  that  shock  tube  end 
wall  boundary  layer  studies  with  the  Wollaston  prism  interferometer 
(Ref.  12),  coupled  with  ultrasonic  and  molecular  beam  techniques  (Ref.  13), 
will  provide  a  firm  experimental  basis  for  transport  properties  of  gases 
in  the  chemical  reaction  region.  Knowledge  of  transport  properties  of 
gases  in  the  temperature  range  1000  <  t  <  8000  K  is  a  continuing  engineer¬ 
ing  need  for  reentry  and  propulsion  systems  design  and  testing. 
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The  immediate  goal  of  the  effort  described  below  has  been  to 
develop  ultrasonic  techniques  for  the  determination  of  thermal  con¬ 
ductivity,  viscosity,  and  diffusion  coefficients  to  support  theoretical 
work  in  the  dense  gas  region.  Early  ultrasonic  pulse  measurements 
showed  that  the  sum  of  thermal  conductivity  and  viscosity  was  useful 
as  a  third  independent  viscosity  measurement  when  a  discrepancy 
between  capillary  and  oscillating  disk  viscosity  measurements  was 
recognized.  The  ultrasonic  pulse  method  was  further  developed  during 
the  past  two  years  and  applied  to  the  measurement  of  transport  prop¬ 
erties  of  argon,  nitrogen  and  air  at  elevated  pressures.  The  sound 
speed  measurements  taken  previously  were  spot  checked  with  a  pulse 
interferometric  technique.  Additional  sound  absorption  measurements 
were  made  in  nitrogen  and  argon  to  fill  in  gaps  in  previous  measure¬ 
ments.  In  addition,  the  frequency  range  of  measurements  with  the 
pulse  technique  was  extended  to  15,  21  and  27  MHz.  New  measure¬ 
ments  of  sound  absorption  and  sound  speed  were  completed  in  air  at 
pressures  up  to  ~500  atm  around  room  temperature.  Transport  and 
thermodynamic  properties  of  air,  nitrogen  and  argon  were  derived 
from  these  measurements  as  described  below.  Further  work  was 
completed  to  show  that  a  technique  based  on  measuring  the  width  of 
acoustic  resonance  lines  in  a  gas  filled  cavity  could  separate  the 
transport  properties  (Refs,  14,  15,  16).  In  addition,  because  the 
measured  quantity  is  frequency,  the  cavity  technique  may  provide 
the  most  accurate  method  for  absolute  determination  of  transport 
properties  known.  The  resonant  cavity  technique  was  set  up  and  mea¬ 
surements  described  below  were  completed  which  prove  the  basic  ideas 
behind  the  measurement.  The  crosshatched  area  on  Fig.  1  denotes  the 
pressure-temperature  region  where  the  cavity  technique  is  applicable. 
Measurements  up  to  1300  K  are  relatively  straightforward.  Extension 
of  the  technique  to  2000  K  with  electrostatic  transducers  seems  feasible 
based  on  our  preliminary  investigations. 

In  the  future,  the  pulse  techniques  (and  the  schlieren  inter¬ 
ferometer)  will  be  used  for  transport  property  measurements  above 
2000  K.  Based  on  our  work  to  date  we  feel  the  accuracy  of  such 
measurements  can  be  kept  to  within  5%.  The  cavity  technic^e  will  be 
used  to  satisfy  transport  property  requirements  below  2000  K  with  a 
precision  of  0.  1%  or  greater.  With  these  three  techniques  we  are  in 
a  position  to  take  data  which  will  extend  the  state-of-the-art  either  in 
the  high  temperature  range  of  3000  <  T  <  13,  000°K,  p  <  200  atm  or 
the  high  density  range  T  <  2000°K  1  <  p  <  1000  atm,  depending  on  the 
most  pressing  needs. 
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SECTION  n 

THEORY  OF  ULTRASONIC  ABSORPTION 

2. 1  PHYSICAL  BASIS  OF  ULTRASONIC  TECHNIQUES 

The  sound  wave  consists  of  small  spatial  variations  in  pressure, 
temperature,  flow  velocity  and  density.  These  gradients  cause  fluxes  of 
heat,  momentum  and  atomic  species  in  opposition  to  the  reversible  fluxes 
in  the  sound  wave.  Gradually  the  energy  of  the  sound  wave  is  lost  in  the 
form  of  heat.  The  losses  eventually  show  up  as  decreases  in  pressure 
amplitude  at  the  receive  transducer.  The  change  in  sound  amplitude 
versus  distance  may  therefore  be  used  to  measure  thermal  conductivity, 
viscosity  and  diffusion  coefficients.  . 

The  sound  wave  is  ordinarily  introduced  by  a  vibrating  solid 
surface  (for  an  example  of  the  opposite  case  see  literature  on  the  Spec- 
traphone  (Refs.  17,  18) ).  The  translational  velocity  of  the  molecules 
adjacent  to  the  solid  surface  moves  in  harmony  with  the  vibrating  solid 
surface  or  transducer.  In  monatomic  gases  the  fluctuations  show  up 
simply  as  an  oscillation  in  pressure,  temperature,  particle  velocity 
(as  opposed  to  translational  velocity),  and  density.  In  gases  with  in¬ 
ternal  modes  of  energy  storage  such  as  vibration,  rotation,  and  elec¬ 
tronic  excitation  the  fluctuations  in  translational  energy  leak  into  the 
internal  modes.  When  the  translational  velocity  falls  below  the  ambient 
value,  the  energy  in  the  internal  modes  begins  to  leak  out.  Since  the 
internal  modes  cause  energy  fluxes  which  oppose  the  reversible  sound 
wave,  absorption  of  sound  energy  results.  Acoustic  losses  due  to  energy 
transfer  between  various  molecular  energy  states  are  referred  to  as  re¬ 
laxation  losses.  The  relaxation  processes  are  of  interest  themselves  in 
transport  theory  for  fast  relaxation  effects  such  as  rotation.  Otherwise 
the  relaxation  effects  must  be  accounted  for  when  one  desires  to  measure 
transport  properties  with  ultrasonic  techniques. 

Intuitively,  one  should  expect  that  a  perturbation  of  the  transla¬ 
tional  velocity  of  gaseous  molecules  travels  with  some  suitable  average 
of  the  thermal  velocity  of  the  molecules.  In  a  monatomic  gas  at  low 
pressures,  the  velocity  of  a  sound  wave  is  ~0.  91  times  slower  than  the 
most  probable  thermal  velocity.  The  sound  speed  measurement,  in  fact, 
is  roughly  equivalent  to  a  gas  thermometer  and  measures  the  translational 
temperature  of  the  gas  it  propagates  through.  One  of  the  most  attractive 
features  of  ultrasonic  techniques  for  transport  property  measurements 
is  that  the  temperature  is  automatically  measured  by  the  sound  velocity. 
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When  the  temperature  is  known  by  other  means  the  sound  speed  and 
its  temperature  or  pressure  variation  are  sufficient  to  determine  all 
the  thermodynamic  properties  of  the  gas.  This  is  very  attractive  at 
elevated  pressures. 

Acoustic  waves  in  a  gas  may  be  free  waves  attenuated  by  con¬ 
version  of  acoustic  energy  to  heat  within  the  gas,  or  wall  waves  which 
propagate  in  a  thin  layer  of  gas  near  a  wall  or  discontinuity.  The  free 
waves  (~1  MHz)  may  easily  be  transmitted  at  distances  up  to  the  order 
of  feet  by  choosing  a  sufficiently  low  acoustic  frequency.  Therefore, 
the  free  waves  are  useful  for  reaching  into  hot  gases  above  the  destruc¬ 
tion  temperature  of  materials  to  determine  transport  properties.  Wall 
waves  decay  very  rapidly  as  they  propagate  away  from  the  wall.  Typical 
values  of  the  attenuation  length  (the  distance  from  the  wall  where  the 
wall  wave  is  attenuated  to  l/e  of  its  original  value)  are  10"^  cm  at  1  kc 
to  10“'^  cm  at  1  MHz.  The  wall  waves  are  the  major  source  of  absorption 
in  cavities  of  the  type  employed  in  the  acoustic  resonance  technique. 

2.  2  EXPRESSIONS  RELATING  TO  TRANSPORT  PROPERTIES, 

SPEED,  AND  ABSORPTION  OF  VARIOUS  SOUND  WAVES 

The  theory  which  relates  the  transport  properties  to  the  measured 
gradients  in  temperature,  pressure,  density  or  concentration  in  any  of 
the  conventional  techniques  for  the  determination  of  transport  properties 
employs  the  linearized  relations  of  fluid  mechanics  and  certain  approxi¬ 
mate  boundary  conditions  (Refs.  19,20,21).  The  linearization  approxi¬ 
mation  is  'more  closely  fulfilled  with  ultrasonic  techniques  than  conventional 
techniques  because  deviations  of  only  1  pt  in  10^  (pulse  techniques)  to  1  pt 
in  10^  (cavity)  of  the  ambient  conditions  are  sufficient  to  conduct  experiments. 

The  reduction  of  the  linearized  equations  of  motion  to  three  wave  . 
equations  can  be  carried  out  with  no  further  approximations  and  is  accom¬ 
plished  in  Appendix  A.  These  three  wave  equations  describe  a  free  wave 
which  corresponds  to  the  sound  wave  employed  in  the  pulse  technique, 
and  two  wall  waves  which  are  responsible  for  the  major  part  of  the  ab¬ 
sorption  in  the  cavity  technique. 

The  free  wave  propagates  with  a  speed  given  to  1  pt  in  10^  by 

c  =  t-zRT  <=■»/«<=  (1) 

and  a  sound  absorption  (in  the  absence  of  relaxation  effects)  equal  to 
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1 

Zp  c 


b> 


4 

3-  n  +  n' 


p 


l/cm 


(2) 


Equation  (1)  serves  as  the  basis  of  temperature  or  thermodynamic  property 
measurements.  Once  the  compressibility  factor  z  is  known  from  other 
measurements  the  specific  heat  ratio  7  or  the  specific  heat  at  constant 
pressure  Cp  may  be  determined  for  sound  speed  and  temperature  deter¬ 
minations.  The  sound  absorption  of  free  waves  may  be  used  to  determine 
any  one  of  the  three  transport  properties,  viscosity  rj ,  bulk  viscosity  rj  ’  , 
or  thermal  conductivity  once  two  of  these  quantities  are  known  from 
independent  measurements.  At  very  high  densities  absorption  of  free 
waves  is  the  only  technique  known  to  the  authors  for  measurement  of 
bulk  viscosity. 

The  velocity  and  absorption  of  wall  waves  depend  on  either  the 
viscosity  or  thermal  conductivity  alone.  The  viscous  wave  is  generated 
when  an  acoustic  field  causes  a  velocity  boundary  layer  near  a  containing 
wall.  The  sound  speed  (c^)  and  absorption  (a  2)  for  this  wave  are  equal  to 


Wimilarly,  the  thermal  wave  which  generated  in  the  vicinity  of  a  wall  due 
to  loww  of  heat  from  the  temperature  fluctuations  in  a  sound  field  has  a 
sound  speed  C2  and  absorption  a  2  given  by  .(neglecting  relaxation  effects 
on  the  specific  heat) 


(5) 

(6) 


The  attenuation  for  both  of  these  wall  waves  is  very  large,  the  wave 
attenuating  by  l/e  in  one  wavelength  which  is  typically  10"^  cm!  The 
viscous  wave  can  be  recognized  as  the  usual  viscoelastic  wave  derived 
in  many  places  for  liquids  (Refs.  22,23).  Basically,  it  is  the  fact  that 
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these  waves  depend  on  only  viscosity  or  only  thermal  conductivity, 
together  with  the  y/1^  frequency  dependence  (as  opposed  to  for 
free  waves)  which  allows  us  to  separate  all  the  transport  properties 
in  the  ultrasonic  cavity  technique. 

The  mathematical  details  of  sound  propagation  are  worked  out 
in  Appendix  A  to  explicitly  demonstrate  that  all  physical  effects  can  be 
rigorously  accounted  for.  This  is  important  to  carry  out  because  the 
cavity  technique  can  apparently  determine  gas  properties  very  accu¬ 
rately.  First,  we  have  shown  in  Appendix  A  that  accurate  usable 
expressions  for  acoustic  losses  in  cavities  may  be  obtained  even  after 
gas  equation  of  state  and  relaxation  effects  have  been  included.  Second, 
we  show  that  the  boundary  value  problem  can  be  solved  accurately 
enough  to  give  0.  1%  transport  property  measurement  and  is  easily 
extended  to  provide  more  accurate  expressions. 


6 


AEDC-TR-71-191 


SECTION  III 

PULSE  TECHNIQUES  AT  HIGH  PRESSURE 

3.  1  EXPERIMENTAL  SETUP  FOR  PULSED  MEASUREMENTS 
AT  elevated  pressures 

Various  fixed  path  techniques  were  developed  under  the  present 
program  in  an  effort  to  obtain  maximum  precision  in  the  sound  absorp¬ 
tion.  The  fixed  path  technique  was  chosen  originally  for  ease  of  operation 
at  elevated  temperature  and  density  (1000°C  <  T,  1000  atm  <  p).  The 
probe  design  and  the  high  pressure  cell  are  consistent  with  operating 
temperatures  at  least  up  to  1000  C.  However,  the  major  effort  has  been 
directed  to  high  precision  sound  absorption  at  temperatures  around  300  K. 
Precisions  between  1  and  0.  1%  have  been  achieved  during  the  past  year. 

The  basic  requirements  for  sound  absorption  and  sound  speed 
measurements  with  the  pulse  technique  are  shown  in  Fig.  2.  A  pulsed 
rf  oscillator  is  used  to  drive  a  transmitter  ultrasonic  transducer.  The 
sound  wave  is  received  by  a  similar  transducer  after  traversing  an 
'o'l  inch  gas  path  and  is  displayed  on  an  oscilloscope.  The  pulse  can 
reflect  from  the  receiver  face,  the  transmitter  face  and  the  resulting 
echo  picked  up  after  traversing  three  times  the  gas  path.  In  practice, 
up  to  4  echoes  may  be  picked  up  before  beam  spreading  or  alignment 
renders  the  pulse  unusable. 

Delayed  oscilloscope  triggers  are  used  to  simultaneously  display 
any  two  echoes  or  an  echo  and  a  sample  of  the  driver  pulse  on  the  dual 
beam  oscilloscope.  An  example  of  a  display  of  the  first  transmitted 
pulse  (the  rounded  trace)  and  the  driver  pulse  on  the  dual  beam  oscil¬ 
loscope  is  shown  in  Fig.  3.  The  delays  are  additionally  used  to  measure 
transit  times.  As  the  pressure  is  changed  a  cycle  of  the  pulse  is  main¬ 
tained  at  a  given  point  on  the  oscilloscope  face  by  varying  the  oscilloscope 
trigger  delay  setting.  The  transit  times  are  measured  this  way  to  the 
nearest  0.  05  |i  sec  in  the  present  work. 

The  frequency  of  the  rf  pulse  is  determined  to  the  nearest  5  cycles 
by  looking  at  the  interference  between  a  cw  signal  fed  into  the  oscilloscope 
in  series  with  the  rf  pulse.  The  cw  oscillator  frequency  is  varied  until 
(Ref.  24)  stationary  parallel  horizontal  lines  appear  across  the  oscilloscope. 
This  means  that  the  cw  and  the  pulse  which  are  incoherent  interfere  at 
constant  phase  points  for  all  times  displayed.  Similar  parallel  lines  will 
appear  when  the  ratio  of  the  pulse  frequency  to  the  cw  frequency  are 
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simple  integers.  This  means  that  high  frequencies  may  be  measured 
by  counting  Some  cw  subharmonics  of  the  pulse  frequency,  A  frequency 
divider  was  used  to  bring  the  1-10  MHz  frequencies  into  the  range  of 
our  1  MHz  counter. 

Finally,  the  temperature  was  measured  with  thermocouples 
located  near  each  transducer.  The  pressure  was  measured  with 
Bourdon  gages,  a  0-1500  psi  Wallace  and  Tiernan  gage  in  the  low 
pressure  range,  and  a  0-7500  psi  Heise  gage  at  elevated  pressures. 

The  temperature  was  determined  to  -O.  03%  and  pressure  to  better 
than  0.  1%. 

3.  2  PULSE  COMPARISON  METHOD 


In  an  effort  to  improve  precision,  measurements  of  the  first 
transmitted  pulse  and  the  driver  pulse  were  used  to  measure  sound 
speed,  density  times  sound  speed,  and  sound  absorption.  The  sound 
speed  was  measured,  as  mentioned  above,  by  keeping  the  position  of 
the  first  pulse  constant  on  the  oscilloscope  face  by  varying  the  oscil¬ 
loscope  trigger  delay.  The  delay  was  digitally  set  to  the  nearest 
0.  1  psec  so  that  the  transit  times  are  accurate  to  the  nearest  _+  0.  05 
|xsec.  Th^  gives  a  precision  of  0.  1%  in  the  worst  case.  There  was 
enough  signal -to -noise  (3000-1)  at  the  elevated  pressures  to  accurately 
locate  the  first  cycle.  However,  at  the  lower  pressures,  ~600  psi,  the 
signal  was  too  low  to  locate  the  first  cycle.  When  a  cycle  was  missed 
it  showed  up  as  a  discontinuity  of  the  transit  time  versus  pressure 
curve.  The  transit  times  had  to  be  corrected  for  phase  shifts  in  the 
receive  and  transmit  transducers.  These  corrections  amount  to  a  half 
period  of  the  rf  frequency;  0.  15  |j.sec  at  3  MHz  and  0.  05  psec  at  9  MHz. 

The  amplitude  information  obtained  with  the  pulse  comparison 
method  may  be  used  to  measure  acoustic  impedance  (pc)  at  low  fre¬ 
quencies  and  free  wave  sound  absorption  (c )  at  elevated  frequencies. 

The  amplitude  of  the  first  transmitted  pulse  is  normalized  by  the  trans¬ 
mitter  cimplitude.  Thus,  fluctuations  in  oscillator  output  and  overall 
gain  are  normalized  out.  The  amplitude  of  the  first  transmitted  pulse 
is  maintained  within  the  range  of  one  of  the  voltage  scales  of  the  oscil¬ 
loscope.  As  the  pressure  increases  less  driver  signal  is  needed  because 
the  acoustic  impedance  of  the  gas  increases  and  the  attenuation  decreases. 
Therefore,  the  ratio  of  transmitted  sound  amplitude  to  driver  pulse 
amplitude  (V^)  increases  as  pressure  increases. 


The  normalized  transmitted  amplitude  (V.)  is  given  by 


Vt  =  A 


\Pt^t) 


f  (p  c,  u)  e 


-cx 
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where  A  is  a  factor  relating  to  the  activity  of  the' transducers,  (pc)  is 
the  acoustic  impedance  of  the  gas,  and  Pt^t  is  l^he  acoustic  impedance 
of  the  transducer.  f(p  c,  <o)  is  a  correction  which  accounts  for  the 
finite  transducer  size  and  bond  thicknesses.  The  function  f(p  c,  co), 
discussed  in  Appendix  B,  is  close  to  unity.  At  3  MHz  the  absorption 
is  small  over  the  pressure  range  of  interest  so  that  is  a  measure 
of  pc.  At  9  MHz  the  absorption  versus  pressure  may  be  determined 
from  Vt  versus  pressure  curves. 

photographs  of  3  MHz  pulses  are  shown  in  Fig.  3  at  the  extremes 
of  pressure  7200  psi  and  200  psi.  It  should  be  noted  that  the  frequency 
and  the  pulse  shape  are  not  changed  by  the  pressure  variation.  A  plot  of 
3  MHz  raw  data  in  nitrogen  is  shown  in  Fig.  4.  The  discontinuities  in 
the  normalized  amplitude  versus  pressure  curves  represent  changes  in 
oscilloscope  scales  used  to  display  the  driver  pulse.  The  relative  gain 
of  the  oscilloscope  preamplifiers  was  calibrated  by  displaying  an  rf 
pulse  of  the  frequency  employed  in  data  taking  on  'two  adjacent  oscil¬ 
loscope  preamplifier  voltage  scales.  The  ratio  of  the  readings  of  the 
same  pulse  on  the  two  voltage  scales  was  taken  as  a  correction  factor 
to  the  data  in  Fig.  3.  This  leaves  only  nonlinearities  in  the  preamplifier 
gain  and  the  transducer  activity  (the  latter  mainly  on  the  transmitter 
side). 


photographs  of  the  driver  pulse  and  received  pulses  at  9  MHz 
for  attenuation  determinations  are  shown  in  Fig.  5.  The  upper  photo¬ 
graph  shows  the  driver  pulse.  The  lithium  niobate  transducers  were 
driven  at  the  third  harmonic  to  obtain  9  MHz.  The  structure  on  the 
driver  pulse  represents  3  MHz  distortion  in  the  Arenberg  rf  oscillator. 
Any  three  MHz  component  on  the  driver  voltage  which  reached  the 
receive  transducer  was  removed  by  the  narrow  bandwidth  (~500  kHz) 
of  the  receive  electronics.  Multiple  echoes  are  also  shown  in  Fig.  5. 
The  multiple  echo  data  were  reduced  to  absorption  to  compare  with  the 
pulse  comparison  values  as  discussed  below.  The  lower  photograph 
shows  the  first  received  pulse  and  the  driver  pulse  used  for  the  pulse 
comparison  attenuation  determinations. 

3.  3  EXPERIMENTAL  RESULTS  WITH  PULSE 
COMPARISON  METHOD 

Measurements  of  p  c  with  3  MHz  ultrasonic  waves  were  made  in 
argon  and  nitrogen.  These  data  have  an  internal  consistency  of  ~0.  1%. 
Additional  measurements  of  sound  absorption  at  9  MHz  were  made  in 
the  same  gases  around  room  temperature.  Sound  speed  determinations 
at  3  and  9  MHz  in  argon  and  nitrogen  were  also  made  simultaneously 
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with  the  amplitude  measurements.  In  principle,  pc,  c  and  absorption 
determinations  allow  both  the  sum  of  viscosity  and  thermal  conductivity 
and  the  thermodynamic  properties  of  the  gas  to  be  determined.  However, 
the  p  c  data  did  not  give  as  accurate  compressibility  values  as  can  be  ob¬ 
tained  from  independent  PVT  measurements  in  the  temperature  ~300  K, 
pressure  (1-500  atm)  range  of  the  present  experiments  (Ref.  25). 

Plots  of  p  c/(p  c)  300  psi  are  shown  for  argon  and  nitrogen  in 
Figs.  6  and  7.  The  data  points  are  consistent  to  about  0.  1%.  However, 
there  is  a  3%  deviation  from  the  "theoretical"  curve  at  the  higher  pres¬ 
sures.  The  theoretical  curve  was  obtained  from  compressibility  data 
(Refs.  26,  27)  and  the  sound  speed  measurement  carried  out  in  the 
present  work.  The  sound  speed  measurements  in  both  gases  are  shown 
in  Fig.  8.  The  accuracy  of  the  sound  speed  measurement  is  better  than 
0.  1%.  Therefore,  the  3%  discrepancy  between  the  experimental  and 
calculated  p  c  at  the  higher  densities  is  probably  due  to  imperfect  re¬ 
flection  from  the  receiver  face  as  discussed  in  Appendix  B.  These 
results  will  be  discussed  further  in  a  paper  in  preparation. 

Sound  absorption  measurements  were  made  at  9  MHz  by  the  pulse 
comparison  technique.  The  results  in  nitrogen  are  plotted  versus  pres¬ 
sure  in  Fig.  9.  Also  shown  are  the  absorption  coefficients  determined 
from  the  amplitude  first  pulse  and  the  first  echo  (top  Fig.  5).  The  first 
pulse-first  echo  data  were  taken  for  completeness.  Absorption  from 
both  types  of  measurement  shows  no  systematic  deviations;  however, 
the  pulse  comparison  data  is  more  precise  in  this  set  of  data. 

The  pulse  comparison  sound  absorption  measurements  require 
the  pc  effects  to  be  divided  out.  The  pc  values  used  were  the  actual 
experimental  results  in  Figs.  6  and  7.  The  theoretical  curve  in  Fig.  9 
is  based  on  other  transport  property  measurements  (Ref.  28)  and  our 
rotational  collision  number  measurements  (Refs.  29,30,  31,32).  All 
the  results  are  consistent  to  ~5%.  The  attenuation  measurements  by 
the  pulse  comparison  technique  have  a  precision  of  1-0.  1%.  The  ques¬ 
tion  of  accuracy  requires  further  comparison  with  other  data. 

Pulse  comparison  measurements  in  argon  are  displayed  in 
Fig.  10.  Here,  the  theoretical  data  is  that  of  Sengers  (Ref.  33).  The 
absorption  measured  in  the  present  work  is  within  a  few  percent  of 
conventional  measurements. 

The  pulse  comparison  technique  seems  capable  of  giving  absorp¬ 
tion  measurements  which  are  self-consistent  to  within  0.  1%  rms  deviation. 
The  authors  feel  that  several  items  in  the  data  reduction  must  be  reexamined 
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before  an  absolute  accuracy  can  be  assigned  to  these  results.  First, 
the  function  f(p  c,  w)  should  be  reexamined.  Second,  the  errors  intro¬ 
duced  through  thermodynamic  properties  employed  must  be  analyzed 
again,  especially  with  regard  to  specific  heats  and  sound  speed. 

Finally,  the  bulk  viscosity  should  be  recalculated  and  uncertainties 
in  rotational  collision  numbers  and  translational  bulk  viscosity  esti¬ 
mated.  The  cavity  measurement  planned  for  the  coming  year  should 
be  useful  for  putting  limits  on  the  bulk  loss  terms.  The  uncertainties 
associated  with  the  above  factors  amount  to  a  few  percent  at  the  present 
time.  Final  calculations  of  these  effects  should  be  completed  in  the 
first  six  months  of  1971. 
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SECTION  IV 

PULSE  ECHO  MEASUREMENTS 


4.  1  HIGH  PRECISION  PULSE  ECHO  FIXED 
PATH  MEASUREMENTS 


Pulse  echo  measurements  were  made  with  several  improvements 
which  came  to  light  during  the  pulse  comparison  work.  Acoustic  non- 
linearities  in  the  gas  near  the  transmitter  and  calibration  of  the  receiver 
electronics  were  accounted  for  experimentally.  This  was  done  by  measur¬ 
ing  attenuation  versus  received  amplitude  and  extrapolating  to  zero  ampli¬ 
tude.  Special  care  was  taken  to  insure  that  the  measurements  were  taken 
at  the  resonant  frequency  of  the  crystal  or  one  of  its  harmonics.  Finally, 
the  frequency  was  measured  with  the  technique  discussed  in  Section  3.  1. 
These  improvements  have  eliminated  most  of  the  scatter  and  excess 
absorption  noted  in  our  previous  work  (Ref.  1). 

Additional  motivation  for  doing  more  pulse  echo  measurements 
came  from  a  desire  to  obtain  a  check  on  the  pulse  comparison  measure¬ 
ments.  The  absorption  from  the  pulse  echo  method  is  absolute  compared 
to  the  relative  values  obtained  from  the  pulse  comparison  technique.  In 
addition,  the  precision  of  the  pulse  echo  measurements  is  about  0,  1%  with 
the  current  improvements.  Finally,  a  technique  based  on  the  interference 
of  the  first  pulse  and  first  echo  was  used  to  obtain  sound  speeds.  This 
technique  eliminates  uncertainties  due  to  finding  the  first  pulse  and  phase 
shifts  in  the'transducers. 

4.  2  EXPERIMENTAL  TECHNIQUE  FOR  PULSE  ECHO 
MEASUREMENTS 


The  experimental  apparatus  for  pulse  echo  measurements  can  be 
seen  by  reference  to  Fig.  2.  A  10  to  1800  cycle  acoustic  pulse  is  initiated 
by  driving  a  lithium  niobate  transducer  with  a  pulsed  rf  oscillator  as  before. 
The  acoustic  wave  travels  through  the  gas  and  is  detected  by  the  receive 
transducer.  However,  most  of  the  energy  (~0,  9999  to  0.  96  in  the  pressure 
range  1-500  atm)  is  reflected  from  the  face  of  the  receive  transducer.  The 
reflected  wave  then  reflects  from  the  face  of  the  transmitter  and  is  detected 
after  traversing  the  path  three  times.  The  echo  is  reduced  in  amplitude  . 
compared  to  the  first  pulse  by  0,  998-0,  921  due  to  reflection  losses  and  by 
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traversing  an  extra  two  path  lengths. 


Measurements  of  frequency,  pressure  and  temperature  are  carried 
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out  as  described  in  the  pulse  comparison  section.  The  absorption  is 
determined  by  measuring  the  amplitude  of  the  first  echo  and  the  first 
pulse  on  a  photograph  of  the  oscilloscope  face.  A  typical  photograph 
is  shown  in  the  lower  part  of  Fig.  11.  The  absorption  is  given  by 


a 


(8) 


where  Aj  and  A^  are  the  amplitudes  of  the  first  pulse  and  the  first  echo 
(corrected  for  reflection  losses),  i  the  probe  separation,  and  a  the 
absorption  coefficient  of  the  ultrasonic  waves.  We  have  adopted  other 
improvements  over  previous  pulse  echo  measurements  which  become 
necessary  as  accuracies  of  1-0.  1%  are  approached.  We  have  mentioned 
gross  heating  effects  due  to  using  too  much  amplitude  in  the  sound  wave 
in  previous  monthlies.  However,  at  lower  amplitudes  there  are  more 
subtle  finite  amplitude  effects  which  simply  lead  to  high  absorption. 

These  effects  must  be  continually  monitored  for  several  reasons.  When 
we  double  the  pressure  leaving  the  frequency  and  driver  amplitude  alone 
the  signal  which  enters  the  gas  is  doubled.  In  addition,  the  average 
amplitude  through  the  gas  path  varies  through  the  gas.  Amplitude  losses 
of  up  to  a  factor  of  100  are  used  to  obtain  accurate  absorption.  Therefore, 
in  the  lower  useful  pressure  range,  where  the  absorption  is  large,  one 
must  be  careful  of  finite  wave  effects  near  the  transmitter.  As  a  result, 
we  now  check  the  attenuation  versus  amplitude  at  constant  frequency  to 
detect  any  possible  finite  wave  effects.  An  example  of  the  type  of  be¬ 
havior  observed  is  shown  in  Fig.  12.  The  extrapolated  values  of  the 
absorption  at  zero  voltage  on  the  transmitter  transducer  is  very  close 
to  the  true  attenuation  as  noted  in  Fig.  12. 

4.  3  EXPERIMENTAL  TECHNIQUE  FOR  PULSE  INTERFERENCE 
SOUND  SPEED  MEASUREMENTS 

At  each  pressure  where  a  pulse  echo  attenuation  was  determined, 
the  sound  speed  was  measured  by  lengthening  the  driver  pulse  width  until 
the  first  pulse  and  the  first  echo  overlap..  The  pulse  and  the  echo  inter¬ 
fere  giving  amplitudes  ranging  between  the  sum  and  the  differences  of 
the  individual  pulse  amplitudes.  In  the  upper  portion  of  Fig.  11  the  first 
pulse,  the  first  echo  and  a  region  of  overlap  with  destructive  interference 
are  shown.  When  the  phase  shift  between  the  first  pulse  and  the  first 
echo  in  the  gas  is  2  tt  n,  the  sum  of  the  two  signals  appears  on  the  os¬ 
cilloscope.  On  the  other  hand,  when  the  phase  shift  across  the  gas  path 
is  2  TT  (n  +  l/2)  the  two  signals  cancel.  The  measurement  of  the  frequency 
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difference  between  two  adjacent  maxima  or  minima  is  a  measure  of  the 
wavelength  in  the  gas.  Since  the  absorption  in  the  high  pressure  cell  is 
large,  maxima  and  minima  due  to  cavity  resonances  do  not  interfere 
with  the  interaction  of  the  two  signals. 

In  practice,  the  frequency  difference  for  10-100  maxima  is  swept 
through,  starting  with  the  rf  oscillator  on  a  frequency  f  j  and  finishing  with' 
frequency  f2-  n  minima  are  counted  off.  The  frequency  shift  between 

successive- minima  at  the  average  frequency  f  =  - - -  is  obtained  by 

dividing  Af  by  n.  The  frequency  shift  varies  linearly  with  average  frequency 
due  to  phase  shifts  in  the  transducer  (Ref.  34)  except  at  the  resonance.  The 
frequency  shift  between  successive  minima  at  the  resonant  frequency  is  ob¬ 
tained  by  fitting  a  straight  line  to  the  frequency  shift  versus  average  fre¬ 
quency  and  picking  a  frequency  shift  (Af  )  at  the  transducer  resonant  fre¬ 
quency  from  the  fit. 

The  sound  speed  is  then  obtained  from  the  transducer  separation  (i  ) 
and  the  measured  frequency  shift  Af^  from 

c  =  2f  Af  .  (9) 

o 

An  example  of  the  first  sound  speed  measurements  in  helium  made  by  this 
technique  is  shown  in  Fig.  13.  Also  shown  is  a  line  based  on  some  inter¬ 
ferometer  measurements  (Ref.  35)  of  sound  speed  in  helium.  The  results 
agree  well  within  0.  1%.  The  errors  of  the  frequency  determination  lead 
to  0.  05%  in  the  individual  sound  speed  measurements.  The  uncertainty  in 
path  length  leads  to  0.  1%  error  in  these  measurements.  The  path  length 
errors  are  systematic  and  do  not  affect  the  relative  accuracy  of  the  sound 
speed  versus  pressure  determination. 

The  interference  method  has  several  advantages  for  sound  speed 
determination  over  the  pulse  overlap  technique  used  to  date.  First,  fre¬ 
quency  of  the  signal  may  be  measured  more  accurately  than  the  transit 
time.  Second,  a  minimum  is  much  easier  to  detect  accurately  than  the 
first  cycle  of  a  pulse.  Finally,  the  problem  of  calibrating  or  synchronizing 
the  triggers  of  the  first  pulse  and  first  echo 'Channels  does  not  arise  in  the' 
interferometric  technique.  Since  the  interference  technique  employs  a 
different  electronic  principle  than  the  usual  transit  time  measurement,  the 
agreement  between  these  two  results  constitutes  a  consistency  test  of  the 
sound  speed  measurements. 
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4.4  RESULTS  OF  PULSE-ECHO  TECHNIQUES 


Sound  speed  was  measured  between  1  and  500  atm  in  air. and 
argon  with  the  pulse  overlap  technique.  The  sound  speed  determined 
by  the  two  techniques  agree  with  each  other  to  within  0.  1%.  The  sound 
speed  versus  pressure  in  air  and  argon  are  plotted  in  Fig.  14.  These 
results  compare  well  with  independent  measurements  also  shown  in 
Fig.  14  (Refs.  36-40).  Accurate  sound  speeds  are  important  because 
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port  properties. 


is  necessary  to  reduce  the  sound  absorption  to  trans- 


The  major  effect  of  thermodynamic  variables  occurs  in  the  thermal 
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conductivity  terms  of  the  sound  absorption  — r: .  The  working  relations 
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for  specific  heat  and  specific  heat  ratio  are,  respectively, 

7  =  (K.J.)  c^ 
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(11) 


where  the  coefficient  of  thermal  expansion  p  and  the  isothermal  com¬ 
pressibility  K-j-  are  obtained  from  compressibility  data  and  c  is  the 
measured  sound  speed.  No  independent  specific  heat  data  are  used. 


The  major  sources  of  error  are  evident  in  Eqs.  (10)  and  (11).  The 
error  in  7  is  twice  the  error  in  sound  speed  (e  j).  Subtracting  from  (1) 
and  squaring  leads  to  a  total  error  of  6  e  j  in  (7-!)^.  The  precision  of  the 


quantities 


C 

P 


therefore  from  Eq. 


(11)  is  4  e  j. 


Since  the  error  in 


compressibility  data  is  about  0.  1  e  j,  the  quantity  which  determines  the 
accuracy  of  the  thermodynamic  inputs  to  the  reduction  of  sound  absorption 
to  transport  properties  is  therefore  the  sound  speed. 


The  sound  speed  measurements  in  air  are  especially  important 
because  of  the  shortage  of  data  on  this  gas.  In  addition,  the  7  and  Cp 
deduced  from  the  sound  speed  and  compressibility  data  for  reduction  of 
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high  pressure  absorption  measurements  in  air  are  plotted  in  Fig.  15. 

These  thermodynamic  data  are  not  available  elsewhere  and  are  of  value 
in  themselves. 

The  sound  absorption  data  in  air  were  used  to  check  the  con¬ 
sistency  of  some  recent  tables  of  air,  viscosity  and  thermal  conductivity 
(Refs.  41,42)  versus  pressure.  This  was  accomplished  by  subtracting 
the  rotational  relaxation  losses,  and  the  translational  bulk  viscosity"^ 
losses  from  the  total  absorption  to  give  the  absorption  due  to  transport 
properties.  Once  this  quantity  is  obtained  a  consistency  test  of  conven¬ 
tional  measurements  can  be  carried  out.  In  the  case  of  air  at  elevated 
densities  the  viscosity  seems  to  be  reliable,  while  the  thermal  conductivity 
measurements  seem  to  have  deviations  up  to  10%.  Therefore,  the  contri¬ 
butions  due  to  viscosity  may  be  subtracted  to  provide  values  of  \.  This 
comparison  is  shown  in  Fig.  16  and  discussed  below. 

Before  discussing  the  comparison  of  the  measurements  carried 
out  under  the  present  program  with  the  correlation  results  of  Refs.  41 
and  42,  the  various  contributions  to  the  sound  absorption  should  be  re¬ 
called.  The  most  interesting  component  of  sound  absorption  in  a  gas 
such  as  air  is  that  due  to  transport  properties  (q^)  given  by 


P  a 


c 


2 


(0 


1 

2 


(.12) 


where  p  is  the  density,  c  sound  speed,  Cp  specific  heat,  y 
ratio,  rj  viscosity  and  \  thermal  conductivity. 


specific  heat 


In  addition,  for  a  polyatomic  gas  such  as  air  there  is  a  contribution 
(about  20%  of  the  total)  due  to  rotational  relaxation. 


^  °  rot  _  ^  Tr_  y-1  _  /  ^  o  \ 
2  "c8  C  i\RT/ 

w  p  '  ' 


(13) 


This  quantity  can  be  calculated  once  the  rotational  relaxation  number 
is  known.  Measurements  of  rotational  relaxation  times  have  been  carried 
out  in  a  muffle  tube  at  Panametrics  and  more  recent  theory  (Refs.  41,42) 


*The  temperature  of  the  present  measurements  is  such  that  vibrational 
relaxation  is  less  than  a  0.  1%  contribution  even  at  elevated  pressures. 
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and  experiment  (Refs.  43,44)  confirm  the  correctness  of  these  rotational 
relaxation  measurements.  The  rotational  collision  numbers  measured  at 
Panametrics  were  used  to  calculate  the  rotational  bulk  viscosity  as  out¬ 
lined  in  a  previous  final  report  (Ref.  1). 

Finally,  at  elevated  densities  there  is  a  translational  relaxation 
time  due  to  multiple  collisions.  The  bulk  viscosity  may  be  estimated 
from  the  relation 

4 

n’  =  I  Xp^  (14) 

m 

where  X  is  the  Enskog  factor,  as  discussed  in  Ref.  1.  Since  the  bulk 
viscosity  is  less  than  20%  of  the  total  absorption  over  most  of  the  pres¬ 
sure  range  of  the  above  experiments,  the  calculated  corrections  to  the 
total  absorption  due  to  the  bulk  viscosity  degrade  the  accuracy  of  the 
measurements  by  at  most  3%. 

The  absorption  values  are  in  fair  agreement  with  the  tables  of 
Refs.  41  and  42  as  can  be  seen  in  Fig.  16,  However,  at  the  higher 
pressure  the  ultrasonic  absorption  measurements  seem  to  be  about  3% 
lower  than  the  correlation  tables.  Assuming  that  this  is  mainly  due  to 
a  discrepancy  in  thermal  conductivity  we  eliminated  the  viscous  losses 
using  the  viscosity  tables  of  Ref.  41.  The  thermal  conductivity  calculated 
from  the  sound  absorption  measurements  is  less  than  the  thermal  con¬ 
ductivity  from  tables  in  the  correlation  by  6%.  This  may  represent  the 
errors  accumulated  in  Ref,  42  by  the  fact  that  Lo,  Carroll  and  Stiel 
were  forced  to  use  a  correlation  based  on  nitrogen  measurements  be¬ 
cause  existing  high  pressure  thermal  conductivity  data  in  air  was  suspect. 

The  sound  absorption  in  argon  is  plotted  in  Fig.  17.  The  con¬ 
sistency  of  this  data  is  about  0.  1%.  In  addition,  the  sound  absorption 
calculated  from  Sengers'  paper  (Ref.  33)  is  in  good  agreement  with  the 
above  measurements.  The  bulk  viscosity  is  less  than  20%  of  the  total 
absorption  and  was  calculated  as  described  above. 
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SEC  TION  V 
CAVITY  TECHNIQUE 

5.  1  GENERAL  DESCRIPTION  OF  MEASUREMENT 

It  is  generally  known  that  sound  absorption  has  contributions 
which  are  proportional  to  viscosity  and  thermal  conductivity.  The 
usefulness  of  ultrasonic  pulse  absorption  data  for  providing  experi¬ 
mental  tests  of  transport  theory  in  high  temperature  (1000-  10,  000°K) 
gases  has  been  amply  demonstrated  (Refs.  1,2,3).  However,  because 
the  effects  of  thermal  conductivity  and  viscosity  could  not  be  separated 
and  because  of  the  difficulty  of  obtaining  accuracies  better  than  a  few 
percent,  the  pulse  technique  is  not  competitive  with  conventional  tech¬ 
niques  in  the  lower  temperature  range  (100-2000  K).  However,  a 
resonant  cavity  technique,  which  has  provided  sound  absorption  to 
within  1%,  was  developed  in  I960  (Refs.  14,  15)  and  has  been  success¬ 
fully  employed  at  the  Physikalisches  Institute,  Stuttgart  (Refs.  16,43) 
since  that  time.  The  present  section  describes  modifications  of  the 
cavity  technique  and  theory  which  potentially  will  provide  extremely 
precise  (<0.  1%)  thermal  conductivity,  viscosity,  and  other  transport 
property  measurements  at  temperatures  below  2000  C. 

The  measured  quantities  in  the  acoustic  cavity  technique  are  the 
resonant  frequencies  of  the  cavity  and  frequency  range  over  which  the 
resonance  has  more  than  half  its  maximum  power  known  as  the 
"linewidth.  "  The  frequencies  at  which  the  resonances  occur  are 
proportional  to  the  dimensions  of  the  cavity  and  the  sound  speed. 
Therefore,  in  a  cavity  of  known  dimensions  the  resonant  frequencies 
provide  the  thermodynamic  properties  of  the  gas.  The  linewidth  on  the 
other  hand  depends  on  the  acoustic  losses.  This  measurement  in  a  cavity 
of  known  dimensions  provides  the  thermal  conductivity,  viscosity,  and 
other  transport  and  relaxation  properties  of  the  gas.  However,  all  of 
these  properties  may  be  separated  because  resonances  at  different 
frequencies  depend  on  the  various  properties  in  independent  ways. 

The  basic  experimental  reasons  for  the  precision  of  the  cavity 
measurement  are  the  conversion  of  all  me'asurements  to  frequency  and 
cavity  radius  determinations.  Using  a  local  Loran-C  station  we  were 
able  to  maintain  a  100  kc  standard  oscillator  to  1  pt  in  10^.  In  principle, 
this  means  that  resonant  frequency  can  be  measured  to  1  pt  in  10®  and 
linewidth  to  1  pt  in  10®. 
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In  order  to  realize  some  of  this  precision  in  practice  the  system 
shown  in  Fig.  18  was  developed  to  measure  the  frequency  of  any  part  of 
the  resonance.  The  procedure  may  be  more  easily  understood  if  the 
resonant  frequency  and  the  linewidth  are  redefined  in  terms  of  phase. 


The  phase  difference  between  the  external  signal  driving  the  cavity  and 
and  the  sound  field  within  the  cavity  increases  by  180°^-^j  as  one  sweeps 


the  frequency  through  a  resonance.  The  resonant  frequency  is  defined 
as  the  90  phase  shift  point.  The  linewidth  is  the  difference  between  the 
frequencies  at  the  45  and  135  phase  points. 


A  phase  lock  amplifier  has  two  inputs,  a  reference,  and  a  signal 
channel.  When  the  reference  and  the  signal  are  in  phase,  the  phase 
lock  amplifier  has  a  maxima.  When  the  reference  and  signal  channels 
are  90  out  of  phase  the  phase  lock  output  is  zero.  First  plots  of  phase 
and  amplitude  are  made  for  the  resonance  being  measured  as  shown  in 
Fig.  25  for  a  typical  cavity  mode.  This  determines  the  resonant  fre¬ 
quency  to  0.  3  Hz  (pt  in  10^)  and  the  location  of  the  45°  frequency  to  1  pt 
in  10^.  Sending  the  external  oscillator  through  the  phase  shifter  on  the 
phase  lock  amplifier  preset  to  shift  the  reference  channel  45°  a  null  is 
obtained  between  the  external  electronic  oscillator  and  field  in  the  cavity 
for  the  45  point  at  a  frequency  f.  A  calibrated  phase  shift  of  90  is 
inserted  between  the  external  oscillator  and  the  reference  channel  of 


the  phase  lock  amplifier.  The  frequency  of  the  external  oscillator  driving 
the  cavity  is  then  increased  until  a  null  is  obtained  at  a  frequency  This 
is  the  135  point  to  within  the  precision  of  the  original  45  phase  point 
frequency  measurement.  However,  because  of  the  line  symmetry,  the 
difference  between  fi  and  fo  is  the  true  linewidth  to  the  order  of 

(  — I  ~  1  pt  in  10^  provided  the  90°  external  phase  shift  is  accurately 


calibrated.  The  use  of  a  phase  lock  amplifier  to  change  the  phase  deter¬ 


minations  into  a  null  allows  the  linewidth  measurement  to  be  carried  out 


with  a  potential  precision  limited  only  by  the  gain  of  the  amplifier  and  the 
harmonic  distortion  of  the  signal.  When  other  experimental  parameters 
such  as  temperature,  pressure  and  frequency  control  and  cavity  dimen¬ 
sion  are  analyzed  for  the  experimental  apparatus  described  below  are 
taken  into  account  a  precision  of  0.  1%  in  the  linewidth  and  0.  001%  in  the 
sound  speed  are  obtained. 


As  discussed  in  the  Introduction  and  Appendix  A,  the  theoretical 
reason  for  the  potential  precision  of  the  cavity  technique  is  the  existence 
of  very  short  range  thermal  and  viscous  waves  whose  losses  are  confined 
to  within  0.  001  inches  of  the  cavity  wall.  Assuming  only  that  the  Navier- 
Stokes  equations  may  be  linearized  (good  to  ~1  pt  in  10^  in  the  present 
case)  the  equations  of  motion  of  the  gas  reduce  to  a  set  of  three  wave 
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equations,  the  usual  free  wave,  a  thermal  wave  and  a  viscous  wave. 
However,  because  of  the  short  range  of  the  thermal  and  viscous  waves 
the  boundary  conditions  on  the  side  wall  and  end  wall  can  be  applied 
separately  to  an  appropriate  approximation.  This  simplification  allows 
very  accurate  expressions  to  be  derived  with  a  minimum  of  labor. 

5,  2  SEPARATION  AND  MEASUREMENT  OF  TRANSPORT 
PROPERTIES 

When  the  cavity  in  Fig.  18  is  driven  at  a  frequency  where  the 
waves  reflected  from  the  cavity  walls  reinforce,  the  signal  on  the  receive 
transducer  goes  through  a  maxima.  The  pressure  distribution  for  the 
second  longitudinal  mode  is  indicated  by  the  density  of  dots  in  Fig.  18. 

In  order  to  excite  the  cavity  into  a  given  mode  the  frequency  of  the  driver 
is  adjusted  to  a  value  close  to  the  resonant  frequency  of  the  desired  mode. 
Since  the  acoustic  wavelengths  are  approximately  the  cavity  dimensions 
in  the  acoustic  case  the  resonant  frequency  is  not  unduly  sensitive  to  cavity 
dimensions. 

A  physical  insight  to  the  separation  of  viscosity  and  thermal  con¬ 
ductivity  can  be  obtained  by  examining  the  particle  motion  in  the  first 
longitudinal  and  radial  modes  which  are  illustrated  in  Figs.  19  and  20, 
respectively.  The  longitudinal  mode  has  particle  motion  parallel  to  the 
axis  of  the  cylinder.  This  leads  to  viscous  losses  only  on  the  cylindrical 
Walls.  The  radial  mode  has  particle  motion  in  the  radial  direction.  Thus, 
viscous  losses  are  restricted  to  the  end  walls.  Because  of  the  difference 
in  the  viscous  losses  between  the  two  modes,  the  thermal  conductivity 
and  viscosity  may  be  obtained  from  measurements  of  two  types  of  modes. 
Similar  conclusions  can  be  drawn  about  compound  modes  where  the 
particle  motion  has  both  longitudinal  and  radial  contributions.  In  addition, 
as  illustrated  in  Figs.  18  and  19,  the  temperature  distribution  is  different 
in  the  various  modes  so  that  the  losses  due  to  thermal  conductivity  also 
may  be  separated.  The  frequency  at  which  the  first  longitudinal  and  first 
radial  (fj  and  respectively)  modes  are  excited  are  given  in  terms  of 
the  length  (d)  and  radius  (r)  of  the  cavity  as 


longitudinal  radial"''  ratio 


*''7  ,  is  the  root  of  the  equation  — -  =  0. 

1  dr 
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The  third  relation  may  be  used  to  eliminate  the  cavity  length  in  favor 
of  the  radius  and  the  measured  frequencies. 


The  separation  of  the  transport  properties  requires  the  mea¬ 
surement  of  three  resonances.  This  is  because  the  linewidth  of  any 
resonance  depends  on  three  independent  quantities,  the  usual  free  space 
absorption  coefficient  6^,  the  longitudinal  wall  losses  6^^,  and  the 
radial  wall  losses  measurement  of  the  first  longitudinal 

^Lil'  second  longitudinal  5  l2»  first  radial  mode  linewidths 
gives  three  equations  in  6^,  ^lvv*  ^Rw*  Alternatively,  the  first  com¬ 
pound  mode  can  replace  one  of  the  above  modes  with  a  linewidth  6 
and  a  wall  loss 


Once  the  linewidths  of  the  first  and  second  longitudinal,  first 
radial  and  first  compound  modes  have  been  measured  we  have  the  fol¬ 
lowing  four  relations 


1st  longitudinal  mode 


2nd  longitudinal 
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Any  pair  of  these  equations  may  be  used  to  eliminate  the  free  wave  loss 
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in  favor  of  the  measured  linewidths  6^1,  ^R1  ^cl' 

measured  frequencies  fj,  f2,  f3  and  f^.  The  bulk  viscosity  so  deter¬ 
mined  should  be  accurate  to  better  than  1%  if  the  linewidth  measure¬ 
ments  are  good  to  0.  01%.  6^  may  then  be  used  to  eliminate  the  free 

wave  losses  from  the  measured  linewidths.  It  is  important  to  note  that 
the  free  wave  losses  were  chosen  to  be  less  than  a  1%  contribution  to  the 
linewidths.  Therefore,  elimination  of  6^  does  not  introduce  significant 
errors  into  wall  losses. 

The  relations  for  the  wall  losses  in  terms  of  Prandtl  number  and 
viscosity  derived  by  Fritsch  are  given  in  Table  1.  These  are  the  first 
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Table  1.  Wall  losses  for  various  mo 
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s  in  the  first  approximation 


R  +  (1  +  R) 


1  +  (1  +  2R) 


7-1 


1/2 


+2  (l-p^)R  +(1-2R) 


1/2 


AEDC-TR-71-191 


approximation  expressions  as  explained  in  Appendix  A.  The  first  ap¬ 
proximation  expressions  are  already  accurate  to  better  than  1%  and 
have  been  demonstrated  experimentally.  In  practice,  the  wall  losses 
would  be  corrected  for  p  c  and  thermal  losses  to  the  cavity  walls,  line 
shape,  and  cross  coupling  of  the  wall  and  free  waves,  Preliminary 
examination  shows  these  effects  to  be  less  than  0.1%  in  magnitude  for 
our  cavity  at  1  atm  and  should  become  smaller  (except  pc  effects)  at 
elevated  pressures. 

As  an  example  of  how  the  transport  properties  are  obtained  we 
will  write  down  the  expressions  for  the  longitudinal  and  compound  mode 
loss 


Dividing  these  two  equations  one  obtains  a  relation  for  the  Prandtl  number 
in  terms  of  the  experimentally  determined  quantities  and 


The  viscosity  may  be  obtained  by  subtracting  these  two  equations,  giving 


Measurements  of  the  first  two  longitudinal  modes,  the  first  radial  mode 
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and  the  first  compound  mode  have  been  demonstrated  in  our  laboratory. 

The  acoustic  cavity  technique  has  the  unique  capability  of  determining 
viscosity,  thermal  conductivity  and  relaxation  effects  separately  in  the 
same  gas  sample.  In  addition,  the  method  offers  potentially  greater 
accuracy  than  existing  techniques. 

5.  3  PRELIMINARY  VERIFICATION  OF  CAVITY  TECHNIQUE 

Measurements  intended  to  proof  the  experimental  approach  to 
gas  cavity  measurements  in  this  laboratory  were  made  in  argon  at  1  atm 
and  ~300  K.  The  cavity  employed  for  these  measurements  is  shown  in 
Fig.  21.  A  0.25"  diameter  cylinder  of  pzt  axially  polarized  material 
with  a  length  corresponding  to  a  1.  5  MHz  thickness  resonator  was  silver 
epoxied  to  a  0.  0635  cm  steel  plate.  The  pzt  transducer  was  backed  with 
a  copper  cylinder.  The  sound  generated  by  the  transducer  reached  a 
similarly  mounted  transducer  by  a  path  through  the  gas,  as  well  as  an 
acoustically  short  circuiting  path  through  the  cavity  walls  (see  Fig.  18). 

The  main  feature  of  this  cavity  was  that  all  the  cavity  walls  were  flush 
stainless  steel.  The  transducers  are  external  to  the  cavity. 

A  diagram  of  the  experimental  apparatus  employed  for  the  cavity 
results  obtained  to  date  is  shown  in  Fig.  18.  A  Loran  receiver  was 
built  and  through  Boston  College  a  100  kc  master  oscillator  which  can 
hold  1  part  in  10®  was  set  up.  The  master  oscillator  frequency  is 
locked  to  the  Loran  signal  frequency  using  a  Lissajous  pattern  on  an 
oscilloscope.  The  master  oscillator  is  used  as  the  clock  in  a  1  MHz 
counter. 

The  frequency  of  the  working  oscillator  which  drives  the  cavity 
is  monitored  with  the  1  MHz  counter.  The  first  step  in  making  measure¬ 
ments  consists  of  sweeping  the  frequency. of  the  working  oscillator  through 
the  frequency  range  of  a  resonance  with  a  vacuTxm  in  the  cavity.  The 
signal  through  the  cavity  under  these  conditions  consists  of  the  short 
circuiting  of  sound  through  the  cavity  walls.  Since  this  is  not  the  desired 
signal,  part  of  the  output  from  the  working  oscillator  is  sent  through  a 
phase  shifter  and  an  attenuator  which  are  set  to  cancel  out  the  short 
circuiting.  The  cavity  is  then  filled  to  the  desired  pressure  of  gas  with 
the  attenuator  and  phase  shift  set  to  cancel  out  short  circuiting  and  the 
frequency  range  of  the  resonance  swept  through  again.  A  trace  of 
amplitude  versus  frequency  such  as  shown  in  Fig.  26  (for  the  second 
longitudinal  mode)  is  obtained.  This  signal  represents  the  resonant 
sound  field  within  the  cavity. 

The  measurement  of  the  second  longitudinal,  first  radial  and  first 
compound  modes  are  completed  in  argon  at  300  K.  In  addition, measurements 
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of  the  first  longitudinal  mode  were  made.  However,  it  was  found  that 
a  membrane  resonance  in  the  end  walls  of  the  cavity  interfered  with  the 
lower  part  of  the  frequency  span  of  the  first  longitudinal  mode.  This  is 
shown  in  Fig.  22,  where  a  plot  of  the  phase  and  amplitude  versus  fre¬ 
quency  is  given.  The  skewness  of  this  graph  shows  the  interference. 

A  more  revealing  plot  is  Fig.  23,  where  the  amplitude  of  the  sound  in 
the  cavity  is  plotted  as  a  function  of  phase  difference  between  the  driver 
and  the  receiver  on  polar  paper.  Ideally,  this  should  be  a  circle  (Ref.  44). 

An  examination  of  Fig.  23  shows  that  the  circle  is  distorted  at 
the  low  phase  angles  (next  to  the  low  frequency  end  of  the  resonance). 
Based  on  the  frequency  symmetry  of  the  resonance  about  the  center 
frequency  one  can  simply  fit  the  ••good”  part  of  the  amplitude  phase 
curve  with  a  circle  and  construct  the  lower  part  of  the  resonance.  This 
constructed  curve  is  shown  in  Fig.  22.  The  linewidth  of  the  constructed 
curve  is  close  to  the  theoretical  linewidth. 

Plots  of  amplitude  versus  phase  are  sensitive  indicators  of 
distortion  due  to  interfering  resonances  in  the  experimental  system. 

Such  resonances  can  arise  from  the  short  circuiting  signal  resonating 
in  the  wall  or  resonances  in  the  response  of  the  transducer.  We  have 
found  membrane  resonances  in  the  cylinder  end  walls  particularly 
difficult  in  this  regard.  A  resonance  in  some  member  of  the  cavity 
which  does  not  couple  to  the  gas  can  be  eliminated  (with  considerable 
labor)  by  careful  vacuum  measurements.  However,  membrane  reso¬ 
nances  Contribute  a  phase  shift  to  the  signal  across  the  gas  path  which 
is  not  present  under  vacuum  conditions. 

Polar  plots  of  the  first  radial,  first  compound  and  second 
longitudinal  modes  are  shown  in  Fig.  24.  These  modes  show  no  in¬ 
terfering  resonances  in  these  measurements.  The  amplitude  and  phase 
shift  versus  frequency  plots  for  these  modes  are  plotted  in  Figs.  25-27. 
The  sound  speeds  can  be  determined  from  the  graphs  to  0.  01%  by  con¬ 
structing  a  horizontal  line  through  the  amplitude  curves  and  bisecting  it 
to  find  the  resonant  frequency.  However,  the  dimensions  of  the  cavity 
used  were  only  held  to  0.  05%.  (A  cavity  has  been  supplied  at  no  charge 
to  us  by  Federal  Products,  a  division  of  Esterline  Corporation,  with 
dimensions  held  to  20  p.in. 

The  linewidths  could  be  determined  from  the  graphs  to  about 
1%.  However,  when  actual  data  is  taken  the  center  frequency  is  ap¬ 
proximately  located  with  graphs  such  as  those  given  and  the  45° 
frequencies  are  determined  from  direct  counter  readings  of  frequency 
at  +  45  .  As  pointed  out  before,  the  symmetry  of  the  phase  shift  versus 


25 


AEDC-TR-71-191 


frequency  relation  means  that  errors  in  linewidth  due  to  errors  in  center 
frequency  are  small.  However,  to  minimize  frequency  drift,  etc.  ,  the 
_+  45  points  should  be  measured  as  rapidly  as^ossible.  An  automatic 
method  for  simultaneously  measuring  the  +  45  points  to  be  described 
later,  has  been  developed. 

The  data  taken  is  summarized  in  Table  2.  Sound  speeds  derived 
from  all  three  modes  agree  with  calculated  values  and  other  measure¬ 
ments  to  0.  1%.  Thermal  conductivity  and  viscosity  values  were  derived 
from  the  second  longitudinal  and  first  compound  modes.  These  are  com¬ 
pared  to  accepted  values  in  Table  2.  The  values  measured  with  the  cavity 
technique  agree  with  the  accepted  values  to  well  within  1%.  No  special 
attempt  was  made  to  control  temperature  (+  0.  5°C)  or  pressure  (+  10  mm) 
in  the  above  work.  However,  the  data  taken  indicates  that  the  technique 
is  basically  sound. 


Table  2.  Transport  properties  measured  with  cavity  technique 


1st  Compound 
Linew  idth 

Hz 

2nd  Longitudinal 
Linewidth 

Hz 

Prandtl 

No. 

Viscosity 

Poise 

198.  5 

185.  5 

0.  660 

2.23 

Michels 
Ref.  33 

- 

- 

0.  669 

2.25 
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SECTION  VI 
SUMMARY 


The  acoustic  pulse  techniques  have  been  thoroughly  explored 
for  transport  and  thermodynamic  property  determination.  The  pulse 
techniques  have  the  unique  capability  of  reaching  into  a  hot  gas 
(2000  <  T  <  15,  000)  and  can  make  valuable  contribution  to  thermo¬ 
physical  properties  at  these  temperatures.  Pulse  techniques  were 
further  developed  during  the  present  program  to  provide  specific 
acoustic  impedance  (pc)  and  attenuation  to  precisions  of  1-0.  1%  from 
amplitude  measurements.  Using  these  techniques  measurements  were 
carried  out  in  argon,  helium,  nitrogen  and  air  at  temperatures  around 
300  C  and  pressures  up  to  500  atm.  Data  reduction  of  air  measure¬ 
ments  was  carried  to  completion.  Thermal  conductivity  and  specific 
heat  ratio  of  air  up  to  ~500  atm  are  reported.  Measurements  in  the 
other  gases  agree  with  attenuation,  acoustic  impedance  and  sound 
speed  calculated  from  independent  data.  Further  analysis  of  this 
data  for  publication  is  in  progress. 

A  new  technique  capable  of  providing  highly  precise  thermal 
conductivity  and  viscosity  determinations  was  developed.  This  tech¬ 
nique  has  been  employed  for  several  years  to  make  1%  sound  absorption 
determinations.  We  have  extended  the  theory  so  that  it  is  suitable  to 
deduce  0.  01%  thermal  conductivity  and  viscosity  values  from  linewidth 
measurements.  In  addition,  we  have  modified  the  cavity  and  transducer 
in  a  way  more  suitable  for  high  pressure  transport  property  determina¬ 
tions.  Initial  Prandtl  number  and  viscosity  values  taken  in  this  apparatus 
show  that  the  technique  is  working  as  designed. 
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APPENDIX  A 

DERIVATION  OF  EXPRESSIONS  DESCRIBING  SOUND 
PROPAGATION  IN  A  CAVITY 

equivalence  of  wave  EQUATIONS  AND  LINEARIZED  NA VIER-STOKES 


The  basic  equations  which  define  the  thermal  conductivity  and  vis¬ 
cosity  are  the  conservation  of  energy,  momentum  and  mass.  In  addition, 
the  equation  of  state  and  an  internal  energy  relation  must  be  specified  to 
completely  determine  the  flow  of  a  fluid  subject  to  given  boundary  condi¬ 
tions.  The  linearized  equations  of  motion  for  a  pure  dense  gas  including 
real  gas  equation  of  state  and  relaxation  of  internal  modes  are 


energy 


9  U 
9t 


+  p  V 


\  V  T  =  0 


(A.  1) 


mome  ntum 


0  V  12  1 

P  Q  +  VP-  nv  v+(t7'  +  jTj)v(Vv)=0  (A.  2) 


mass 


8p  ,  _  — 

9t  ^o 


0 


(A.  3) 


equation  of  state  p  =  zp  RT 


(A.  4) 


...  au 

internal  energy  -rr" 
dt 


(A.  5) 


Here  the  internal  energy  U  is  defined  in  terms  of  a  complex  specific  heat 
Cp  and  the  translational  temperature  T.  One  could  define  the  temperature 
statistically  (Ref.  A.  1)  and  have  relaxation  effects  show  up  in  the  bulk  vis¬ 
cosity  rj'  or  use  the  nonequilibrium  thermodynamics  formalism  (Ref.  A.  2). 
However,  we  chose  to  define  a  complex  specific  heat  for  brevity  which  leads 
to  identical  answers.  Methods  for  obtaining  complex  specific  heats  are 
worked  out  elsewhere  (Refs.  A.  3,  A.  4).  We  will  simply  carry  through 
the  calculation  to  show  how  relaxation  effects  enter  the  absorption  and 
velocity  relations. 

Equations  (A,  1-A.  5)  constitute  5  equations  in  the  5  unknowns: 
internal  energy  U,  pressure  p,  velocity 'v,  translational  temperature  T, 
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and  density  p.  With  appropriate  boundary  conditions  expressions  for  the 
flow  field  in  terms  of  the  thermal  conductivity  \  ,  viscosity  rj ,  volume 
viscosity*  jr] '  ,  and  specific  heat  can  be  derived.  The  most  common  way 
to  derive  the  desired  expressions  for  attenuation  and  sound  speed  is  to 

assume  solutions  of  the  type  e  '  exist  and  by  substitution  into 

(A.  1-A.  5)  derive  values  for  k  and  a  (Ref.  A.  4).  However,  Fritsche  (Ref.  5) 
was  able  to  reduce  the  above  five  equations  to  a  set  of  four  wave  equations 
with  no  further  approximations.  Proceeding  in  this  way  allows  all  the 
physical  phenomena  to  be  evaluated  easily.  In  addition,  boundary  condi¬ 
tions  may  be  applied  to  the  resulting  wave  equations  in  a  standard  way. 

The  velocity  may  be  separated  into  a  rotational  and  irrotational 
part  in  the  usual  manner 


V  =  V^+VxH^  (A.  6) 

where  (ft  is  the  scalar  potential  and  the  vector  velocity  potential.  The 
momentum  equation  (A.  2)  then  becomes 


/8H 

+  Vx^— 


0 


with  the  new  scalar  potential  §  defined  as 


§ 


9(i) 

9t 


(ft  + 


_E_ 


(A.  7) 


(A.  8) 


The  wave  equation  for  viscous  waves  which  propagate  near  a  wall 
separates  immediately  from  (A.  7)  giving 

=  0  (A.  9) 

The  wave  number  for  these  waves  is 


(A.  10) 


*The  volume  viscosity  in  the  present  formalism  includes  only  many  body 
contributions. 
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The  velocity  pressure  and  density  may  be  eliminated  from  the 
scalar  part  of  the  momentum  and  energy  equation  (A.  8)  using  the  equation 
of  state 


(A.  11) 


(A.  12) 


(A.  13) 


the  internal  energy  equation  (A.  5),  the  conservation  of  mass  (A.  3)  and  the 
definition  of  the  velocity  potential  (A.  6).  The  resulting  equations  are  with 
no  approximations 


V  §  =  0 


(A.  14) 
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7 


o 

7 


iw(-|  n  +  n') 


Ty2  2 
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=  0  (A.  15) 
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(A.  16) 


a  T 

where  9  is  j  ,  and  the  thermodynamic  identities 


P  7 

have  been  used. 


and 
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Each  of  the  equations  (A.  11),  (A.  12)  and  {A.  13)  resemble  wave 
equations;  therefore,  it  makes  sense  to  see  if  the  variables  $,  0  and  6 
can  be  replaced  to  three  quantities  (f  =  0,  1,  2)  which  obey  wave 
equations 

+  k  H.  =0  (A.  18) 

ill 


A  simple  linear  expansion  in  the 

3 

0=5 

3 

0  =  S 
0 

3 

§  =  S 
0 


is  the  most  obvious  route  to  try 

(A.  19) 

a.  H  (A.  20) 

I  I 

b^  (A.  21) 


If  values  of  ,  a^g  and  b^^  can  be  found  then  the  set  of  wave  equations 
(A.  14)  may  be  used  to  apply  the  boundary  conditions  of  the  particular 
problem  at  hand. 

Substitution  of  (A.  19-A.  21)  in  (A.  14-A.  16)  provides  the  desired 
quantities.  Substitution  in  (A.  11)  shows  that  if  bjg  0  k^  =  0  and  the  first 
wave  equation  is  immediately 

=  0  and  k  =  0  (A.  22) 

0  o 


The  remaining  cases  are  for  bjj  =  0  k^  ^0.  After  substitution  of  (A.  19) 
through  (A.  21)  into  (A.  15)  and  (A.  16)  one  obtains  for  a^ 


a 


(A.  23) 


and  fnr  W  the  nnadra^tic  eaiiation 


2 

c 

o 

2 

Cl) 


+ 


4 

■3  TJ  +  T?' 
i(o  p 


P 


C 

o  V 


+  1=0 


(A.  24) 
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Finally,  after  much  laborious  algebra  and  approximation,  the  two 
remaining  wave  numbers  follow  from  (A.  24) 


(0 


1  -  i 


Cl) 


(I" 


P  P  / 


(A.  25) 


_  p  C 

,2  .  —  V 

^2  =  — 


_  ._<o_  U  ^ _ ^ 

^2  y  7-1 


(A.  26) 


The  basic  equations  of  motion  (A.  1-A.  5)  have  been  reduced  to 
three  wave  equations  with  wave  numbers  given  by  (A.  9,  A.  25)  and  (A.  26). 
These  wave  numbers  are  explicit  functions  of  transport  properties,  re¬ 
laxation  times,  thermodynamic  quantities  and  frequency.  The  approximations 
which  lead  to  (A.  25)  and  (A.  26)  from  (A.  24)  are  not  necessary  but  lead  to 
errors  of  only  one  part  in  IQI^  in  kj  and  k2-  In  addition  to  one  part  in  10^, 
the  second  terms  in  (A.  25)  and  (A.  26)  may  be  neglected.  These  terms  repre¬ 
sent  the  coupling  between  wall  waves  and  free  waves  that  one  might  expect 
physically. 

The  waves,  Hj,  H2  and  H3,  related  to  velocity  pressure  and  tem¬ 
perature  through  the  potentials  defined  in  Eqs,  (A.  6-A.  8)  and  the  relations 
(A.  19-A.  21).  Expressions  for  v(Hi,  H2,  H3),  pCHj,  H2)  and  TCHj,  H^) 
can  be  obtained  from  the  relations  (A.  6-A.  8)  and  (A.  17-A.  19),  in  terms 
of  the  known  wave  numbers  kj,  k2,  k^  and  the  quantities  a^  and  a2.  These 
relations  transform  the  boundary  condition  on  pressure,  temperature  and 
flow  velocity  into  conditions  on  the  solutions  Hj,  H2  and  Hj. 

Following  the  procedure  for  separation  of  variables  one  might 
expect  the  first  wave  equation 

VHi  +  k^^  =  0  (A.  27) 


to  have  solutions  of  the  sort 

00 


p=  0 


r)  e 


z 


(A.  28) 


where  y 
and  are 


and  k_^^^  are  to  be  determined  from  the  boundary  conditions 


required  ^  satisfy 


(1)' 


+  k 


(1)‘ 


=  k 


(2) 


(A.  29) 
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One  might  expect  to  proceed  in  a  similar  fashion  with  the  separation 
of  the  wave  equations  for  H2  and  H2.  However,  in  general,  in  acoustic 
resonators,  because  Hj,  H2  and  H3  are  coupled,  the  solutions  of  separated 
wave  equations  cannot  be  made  to  fulfill  all  the  boundary  conditions  on  the 
end  walls  and  side  walls  simultaneously.  Physical  insight  provides  a  unique 
separation  of  the  wave  equations  which  allows  the  boundary  to  be  approxi¬ 
mately  satisfied.  For  instance,  the  separation  for  H3  may  be  made  in  the 
form 


o  p 


J  (7  ^^^r)e  P 
o  p 


(A.  30) 


where  the  first  term  represents  a  thermal  wall  wave  along  the  side  wall  of 

ik_(l)^z 

the  resonator  excited  by  the  Hj  wave  e  .  And  the  second  term 

represents  a  thermal  wave  on  the  end  wall  of  the  resonator  driven  by  the 
free  wave  Jq  (7p^^^r).  A  similar  separation  can  be  effected  for  the  viscous 
waves  (H^). 


In  order  for  these  waves  to  be  solutions  to  the  wave  equations  for 


and 


P  P  2 


(A.  31) 


k  +  7  -  k  ^ 

V  ^1  2 


(A.  32) 


must  hold. 


Fritsch  showed  that  k  k,  and’ 7  ~  7  (k,  or  7  are  the  wave 

P  1  P  P  ^  1 

parameters  for  the  zero  viscosity  thermal  conductivity  case).  This  leads  to 

(2)^  2  (3)  (3) 

7  '  '  ~  k_  and  similar  results  for  7  '  '  and  k  '  Therefore,  the  first 

P  2  P  P 


term  in  Eq.  (A.  30)  is  small  at  any  distance  from  the  side  wall  of  the  cavity 


(—  at  r  =  0.  001).  Similarly,  the  second  term  is  small  at  any  distance  from 
0 


the  end  wall.  This  is  because  the  wall  waves  fall  off  rapidly  as  we  leave  the 
wall.  An  excellent  approximation  may  therefore  be  obtained  by  applying  the 
boundary  conditions  on  the  side  walls  of  the  cavity  neglecting  the  second 
terms  in  (A.  30).  This  gives  traveling  waves  in  both  directions  including 
the  viscous  and  thermal  waves.  A  reflection  coefficient  for  these  waves 


may  be  obtained  at  the  end  walls  using  only  the  second  terms  in  Eq.  (A.  30), 
which  satisfy  the  end  wall  boundary  conditions. 
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A  complete  solution  can  be  built  up  in  the  above  fashion  leading 
to  the  expressions  for  linewidth  and  resonant  frequency  given  in  the 
text.  These  expressions  neglect  stress  and  temperature  fluctuations 
in  the  wall  and  imperfect  reflection  from  the  end  walls.  In  addition, 
some  algebraic  approximations  are  made  in  the  linewidth  expressions. 
However,  the  first  order  solution  is  already  good  to  0.  1%  over  the 
range  of  operating  conditions  anticipated,  200  <  T  <  1000°C,  p  <  500 
atm.  Corrections  to  the  linewidth  can  be  derived  by  a  simple  perturba¬ 
tion  method  to  correct  the  measured  linewidth  to  within  0,  01%  of  the 
first  order  expressions.  These  corrections  will  be  presented  in  a 
later  report. 
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APPENDIX  B 

TRANSDUCER  ANALYSIS  FOR  PULSE  TECHNIQUE 


The  transducers  employed  in  pulse  sound  speed  and  attenuation 
measurements  are  the  key  to  obtainable  precision  and  temperature  range. 
The  transducer  needs  to  be  sufficiently  damped  to  provide  pulses  which 
do  not  interfere  with  each  other.  In  addition,  the  pulses  must  reach 
their  steady  state  value  in  a  reasonable  time  (^ZOpsec),  The  range  of 
frequencies  at  which  the  transducer  has  usable  output  in  the  neighborhood 
of  a  resonant  frequency  increases  as  the  damping  increases.  However, 
output  decreases  with  more  damping.  The  ability  to  resonate  at  harmonics 
also  decreases  as  damping  is  increased.  The  transducer  mounting  is  de¬ 
signed  as  a  compromise  between  these  effects.  In  addition,  the  transducer, 
bonding  and  backup  materials  must  be  consistent  with  the  highest  desired' 
working  temperature.  Finally,  alignment  and  parallelism  of  all  parts  are 
important  if  accurate  estimation  of  phase  shifts  and  interference  in  the 
transducer  body  are  required. 

The  transducer  employed  in  all  the  measurements  reported  is 
shown  in  detail  in  Fig.  IB.  The  acoustically  active  elements  shown  in 
the  lower  part  of  Fig.  IB  are  the  most  important.  A  lithium  niobate 
piezoelectric  crystal  is  bonded  with  silver  epoxy  to  a  zinc  backup.  The 
zinc  provides  a  siifficient  sound  absorption  and  pc  match  to  the  lithium 
niobate  transducer  to  provide  adequate  damping.  The  rise  time  of  the 
lithium  niobate  transducers  with  zinc  backup  plates  is  about  30  cycles  at 
their  3  MHz  fundamental  frequency. 

The  active  transducer  elements  are  bonded  into  the  supporting 
stainless  steel  transducer  body  with  high  temperature  cement  as  shown 
in  the  upper  part  of  Fig.  IB.  The  front  electrode  is  made  by  coating  the 
transducer  assembly  with  silver  epoxy  and  grinding  the  layer  to  about 
0.  003"  thick.  The  transducer  body  slips  into  a  mount  which  maintains 
alignment  and  makes  electrical  contact. 

The  waves  reflected  from  the  interfaces  of  the  various  layers  in 
the  transducer  interfere  causing  delays  and  apparent  attenuation  of  the 
sound  waves.  For  instance,  the  first  pulse  entering  the  transducer  is 
delayed  by  the  transit  time  in  the  front  silver  layer.  In  addition,  if  the 
first  half  cycle  is  not  used  to  determine  the  arrival  time,  interference  of 
the  waves  reflected  from  the  transducer  interface  causes  an  additional 
delay  of  up  to  l/ZX.  Finally,  waves  which  have  entered  the  transducer 
can  come  out  again  and  interfere  with  the  signal  in  the  gas.  Each  reflec¬ 
tion  leads  to  addition  of  this  "reradiated"  signal  and  leads  to  a  modification 
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of  the  signal  in  the  gas.  When  the  transducer  is  operated  at  resonance 
the  reradiated  signal  is  in  phase  with  the  echo.  However,  when  the 
transducer  is  operated  off  resonance  the  attenuation  may  appear  to 
increase  with  the  magnitude  of  the  difference  between  the  operating 
frequency  and  the  resonant  frequency. 


Finally,  when  using  a  transducer  such  as  shown  in  Fig.  IB  for 
pulse  comparison  determinations  of  p  c  and  attenuation  the  correct 
reflection  and  transmission  coefficients  must  be  used.  For  instance, 
the  relation  between  the  sound  wave  incident  on  the  receive  transducer 
and  the  sound  wave  in  the  transducer  involves  the  sum  of  signals  which 
are  multiply  reflected  within  Ae  front  silver  epoxy  layer.  To  a  lesser 
extent,  the  signal  which  leaves  the  transducer  enters  the  backup  layer 
and  the  bond,  subsequently  reflecting  into  the  transducer,  also  contri¬ 
butes  to  the  transducer  response. 


We  will  calculate  the  transmission  coefficient  of  the  front  layer 
of  the  receive  transducer.  The  signal  enters  the  layer  from  the  gas 
with  a  reduction  of  amplitude  T^^  given  by 


P  c 
o  o 


Pl^l  +Po‘'o 


Pl‘^1 


(B.  1) 


where  P  ^  is  the  acoustic  impedance  of  the  gas  and  P 
impedance  of  the  layer. 


is  the  acoustic 


The  wave  in  the  layer  then  reflects  from  the  transducer-front  layer 
interface.  Some  of  the  incident  energy,  namely 


=  T 


ikjf 


ol  12 


(B.2) 


where  i|;  is  the  incident  signal  in  the  gas,  i|;  the  signal  in  the  transducer 
O  I 

and  T  is  the  transmission  coefficient  at  the  transducer-front  layer  inter¬ 
face 


^l"l 


(B.3) 


is  the  wave  number 


in  the  front  layer. 


/ 

/ 
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Considering  the  multiple  reflections  of  a  signal  in  the  front  layer 
as  leaking  some  signal  into  the  transducer  at  each  reflection  we  can  add 
up  the  total  signals  entering  the  transducer.  For  a  train  of  n  pulses  this 
is 


't  +0  ^ol  ’'12 


l+Rj^e 


2  i  nk  f 

ikji  l-(Ri2Rio)"e 


^  "  ^12^10® 


2ikji 


(B.4) 


where  the  reflection  coefficients  at  the  front  layer-transducer  interface 
and  the  front  layer-gas  interface  are 


.  Pl‘=l-P2'=2 

12  -  PjCj  +P2C2 


and 


10 


Pl®l-Po®o 

Po®o 


1. 


(B.5) 


(B.  6) 


The  important  thing  to  note  is  that  the  signal  in  the  transducer  depends  on 
the  frequencies  {via  k^  =  ~^)  and  the  thickness  of  various  layers.  In 

addition,  the  dependence  on  the  acoustic  impedance  of  the  gas  is  much 

P  0^0 

more  complex  than  the  simple  -  dependence  used  to  explain  the 

^  1®1 

principle  of  the  measurement  in  the  text. 

Fortunately,  the  magnitude  of  approximately  1  for  a  gas  and 

Ri2  ~  0.  5.  Therefore,  the  complex  terms  added  to  1  in  Eq.  (B.  4)  are 
small.  However,  one  has  to  be  careful  because  when  the  front  layer  is 
~  l/4\  thick  and  the  acoustic  impedance  of  the  layer  is  the  geometric  mean 
of  the  acoustic  impedance  of  the  transducer  and  the  gas  the  signal  in  the  gas 
is  entirely  transmitted  into  the  transducer!' 

The  important  transducer  parameters  are  listed  in  Table  IB.  The 
thickness  of  the  various  elements  of  the  transducers  are  given.  A  short 

calculation  shows  that  the  front  layers  are  not  an  even  multiple  of  for 

either  of  the  transducers  at  the  operating  frequencies.  In  addition,  the 
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Table  IB.  Transducer  mounting  and  parameters 


Transducer  Parameters 

Thickness  (cm) 

Receive 

Transducer 

Transmit 

Transducer 

Front  silver 
layer 

0.  0071 

0.  0081 

Lithium  niobate 

0.  106 

0.  106 

Bond 

0.  005 

0.  005 

Zinc 

0.  635 

0.  63  5 

Probe  separation  2.4809  +  0.0019  cm 

Material  Parameters  (300°K) 

LiNbO^ 

Zinc 

Silver  Epoxy 

1  atm  10  100  500 

Sound  speed  cm/p,  sec 

(longitudinal 

0.  652 

0.421 

0.  1909 

0.  0353  0.354  3.  80 

X  10-3 

3  MHz 

Wavelength 

9  MHz 

0.  217 

0.  0724 

0.  140 

0.  0463 

0.  0636 

0.  0212 

Acoustic  gm/p-sec  cm^ 

impedance 

3.  025 

3.  00 

0.  5011 

0.0402  0.404  4.31  22.6 

X  10"^ 
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acoustic  impedance  of  the  front  layers  never  reaches  the  geometric  mean 
of  the  acoustic  impedance  of  the  gas  and  the  transducer  in  the  pressure 
range  of  these  experiments.  In  fact,  for  nitrogen 


'p_c_p  c 
2  2  o  o 


<  0.261 


A  complete  description  of  transducer  corrections  is  beyond  the 
scope  of  this  report.  However,  the  information  in  Table  IB  defines  all 
the  necessary  parameters.  Transducer  characteristics  were  also 
measured  using  the  frequency  behavior  of  the  transducer.  These  results 
will  be  reported  more  completely  in  later  papers  on  the  pulse  techniques. 
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USEFUL  TECHNIQUES  IN  THE 
HIGH  TEMPERATURE  REGION 

ELECTRICAL  CONDUCTIVITY 
CONSTRICTED  ARC  PROBES 
ULTRASONIC  ATTENUATION 


PLASMA 


Z  \/Xf 

„ _ 1.2 

HEAT  FIELD 
^  SCHLIEREN 
>APPLICABLE^ 


USEFUL  TECHNIQUES  IN  THE 
DISSOCIATION  REGION _ 

ULTRASONIC  ATTENUATION 

END  WALL  BOUNDARY  LAYER -WOLLASTON 

PRISM  INTERFEROMETER 


DENSITY  EFFECTS 
\  IMPORTANT  v\ 


>o 


NEW  HIGH  PRESSURE 
EXPERIMENTAL  / 


ULTRASONIC  ATTENUATION 
HEAT- FIELD  SCHLIEREN 
END  WALL  HEAT  TRANSFER 
ULTRASONIC  CAVITY 

^  LITHIUM  NIOBATE 
S  ELECTROSTATIC 


600  800 

PRESSURE  ATM 


1000  1200 


Fig.  1.  Useful  pressure-temperature  ranges  of  various 
techniques  of  transport  property  determination 
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Fig,  2,  Schematic  of  experimental  apparatus  for  pulse  technique 
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3  MHz  Driver  Pulse  and  Transmitted 
Pulse  p  =  7200  psi  =  8.  1  Volts 


3  MHz  Driver  Pulse  and  Transmitted 
Pulse  p  =  200  psi  =  420  Volts 


Fig.  3.  3  MHz  relative  amplitude  acoustic  impedance  traces 
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Fig.  4.  Raw  acoustic  impedance  data  in  nitrogen 
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9  MHz  driver  pulse  and  multiple  echo 


9  MHz  transmitted  pulse  and  driver  pulse 


Fig,  5. 


9  MHz  relative  amplitude  acoustic  attenuation  traces 
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NORMALIZED  ACOUSTIC  IMPEDANCE  Ln(pc/p  3QQC3QQ) 


PRESSURE  ATM 


Fig.  7.  Acoustic  impedance  relative  to  300  psi  in  nitrogen 
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SOUND  ABSORPTION  cm 
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Sound  speed  by  interference  of  1st  pulse  and  1st  echo 


l£t  echo 


1st  pulse 


Attenuation  trace 


Fig.  11.  9  MHz  through  transmission  multiple  echo  and  pulse-echo 
interference  traces. 
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SPECIFIC  HEAT  JOULES/gm"K 


SPECIFIC  HEAT  RATIO  y 
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300 


400 


500 


PRESSURE  ATM 

Fig.  16,  Thermal  conductivity  of  air  versus  pressure 


DENSITY  TIMES  TOTAL  ABSORPTION  po(gm/cm 


PRESSURE  ATM 


Fig,  17.  Sound  absorption  versus  pressure  for  argon  by  the  pulse  echo  technique 


SET  TO  NULL  SHORT  CIRCUITING  UNDER  VACUUM 


Fig.  18.  Diagram  of  apparatus  for  acoustic  cavity  technique 
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SCALE  x2 


Fig.  21.  Details  of  acoustic  cavity  used  in  initial  transport 
property  measurement  with  the  cavity  technique 
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Fig.  Z4.  Polar  plot  of  amplitude  versus  phase  lor  first  radial  mode;  second  longitudinal 


mode  and  first  compound  mode. 
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A.  ASSEMBLY 


HIGH  TEMPERATURE 
CEMENT 


V/A  =  STAINLESS  STEEL 


=  SILVER  EPOXY 


Fig.  I B.  Details  of  transducer  assembly  for  pulse  techniques 
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13.  ABSTRACT  Transport  property  measurements  with  ultrasonic  techniques  are 
discussed  over! the  entire  temperature  ( 100-15, 000°K)  and  pressure  (1- 
1000  atm)  ranges.  A  chart  of  various  techniques  for  the  determination  of 
transport  properties  in  the  various  temperature  and  pressure  ranges  is 
given.  Sound  absorption  and  sound  speed  determinations  were  used  to  ob¬ 
tain  thermophysical  properties  of  argon,  helium,  nitrogen  and  air  at  room 
temperature  and  pressures  between  1-500  atm.  Several  variations  of 
pulsed  ultrasonic  techniques  were  used  to  obtain  maximum  precision  in  the 
gas  property  measurements.  Sound  speeds  accurate  to  approximately  0.1  % 
were  obtained  in  all  the  gases.  Thermal  conductivity  and  specific  heat 
ratios  were  deduced  for  air.  A  cavity  technique  was  developed  for  high 
accuracy  (<0.1  %)  viscosity,  thermal  conductivity  and  bulk  viscosity 
determinations.  Initial  tests  were  run  to  proof  the  technique  at  room 
temperature  and  pressure.  Sound  speed,  Prandtl  number  and  viscosity 
measurements  were  made  which  agree  with  accepted  values  to  within  1  %. 
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